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Induced interactions in a superfluid Bose-Fermi mixture 
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We analyse a Bose-Einstein condensate (BEC) mixed with a superfluid two-component Fermi gas 
in the whole BCS-BEC cross-over. Using a quasiparticle random phase approximation combined 
with Beliaev theory to describe the Fermi superfluid and the BEC respectively, we show that the 
single particle and collective excitations of the Fermi gas give rise to an induced interaction between 
the bosons, which varies strongly with momentum and frequency. It diverges at the sound mode of 
the Fermi superfluid, resulting in a sharp avoided crossing feature and a corresponding sign change 
of the interaction energy shift in the excitation spectrum of the BEC. In addition, the excitation 
of quasiparticles in the Fermi superfluid leads to damping of the excitations in the BEC. Besides 
studying induced interactions themselves, these prominent effects can be used to systematically 
probe the strongly interacting Fermi gas. 


The interplay between induced interactions and su¬ 
perfluidity plays an important role in low temperature 
physics. In metals, the phonon mediated interaction be¬ 
tween electrons leads to the formation of Cooper pairs [I], 
and induced electron-hole excitations significantly sup¬ 
press the critical temperature of a BCS superconduc¬ 
tor mm- a prominent theory for high temperature su¬ 
perconductivity is that it is caused by spin fluctuations 
leading to an attractive interaction [1] , and induced inter¬ 
actions are important for understanding the properties of 
liquid helium mixtures [5j. The systems where induced 
interactions are significant often consist of fermionic and 
bosonic degrees of freedom. In cold atom gases, Bose- 
Fermi mixtures have been realised experimentally for 
sympathetic cooling [SHE], molecule formation ora, and 
for studying few-body physics m- The theoretical stud¬ 
ies have focused on mixtures where the Fermi gas is in the 
normal state [l4l - (23j . Recently, an experimental break¬ 
through was reported with the realisation of a mixture of 
superfluid 7 Li and 6 Li gases [23]. This opens up the excit¬ 
ing possibility to experimentally study for the first time 
the role of induced interactions in a Bose-Fermi mixture, 
where both components are superfluid. 

Here, we study a BEC mixed with a two-component 
superfluid Fermi gas in the whole BCS-BEC crossover at 
zero temperature. Using a quasiparticle random-phase 
approximation (QRPA) to describe the excitations in 
the Fermi gas, combined with Beliaev theory for the 
bosons, we show how the fermions give rise to an induced 
frequency/momentum dependent Bose-Bose interaction, 
which diverges at the sound mode of the Fermi gas. This 
results in two qualitatively new effects. First, the dis¬ 
persion relation of the bosons in the BEC is significantly 
changed at frequencies/momenta close to the sound mode 
of the Fermi gas. Second, bosonic excitations are damped 
due to dissipation, as they can excite quasiparticles in 
the superfluid Fermi gas [25]. These effects can be used 
to systematically probe the single particle and collective 


properties of the strongly correlated Fermi gas. 

We consider a gas of bosons with mass m B mixed with 
a two-component (<j J.) gas of fermions with mass mp. 

The populations of the two fermionic states are taken to 
be the same. The Hamiltonian of the Bose-Fermi mixture 
is H = Hb + Hp + Hbf, where 

= X efcQ k«k + ^- ^ F B (g)a ] J c+q aj i; _ q a k m k (1) 

k k,k',q 

is the Bose hamiltonian with e*, = fc 2 /2?n B , 
fc 2 1 . 

HF = J2i 2 m F C k- Ck - + v X M9)4+ qt c k'- q +c k 'iC k t 

kcr k,k',q 

( 2 ) 

is the Fermi Hamiltonian, and 

Hbf = v X ^BF(<7)c k+qo .c k(T a k/ _ q a k ' (3) 

k,k ; ,qcr 

is the Bose-Fermi interaction. The operators a k (c k(T ) 
remove a boson (spin a fermion) with momentum k, V 
is the volume of the system, and we work in units where 
Tl = fc B = 1. In the following, we replace the interactions 
with the corresponding low energy scattering matrices: 
Vb(<7) —> Tb = 47t a B /m B , Vp{q) —> 7f = 47raF/mF, and 
Vbf(<7) —■► 7bf = 27ra B F/m r , where a B , a B F, and of is 
the Bose-Bose, Bose-Fermi, and Fermi-Fermi scattering 
length respectively, and m T = m B niF/(™ B + Trip) is the 
reduced mass. As usual, this corresponds to summing all 
ladder diagrams in a vacuum. 

The presence of the Fermi gas induces an effective in¬ 
teraction between the bosons, since one boson tends to 
attract/repel fermions giving rise to a local change in the 
fermion density, which is felt by the second boson. Com¬ 
bined with the direct Bose-Bose interaction, this results 
in total interaction 

v(q,uj) = T B + Tg F x(q,u). 


(4) 
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FIG. 1. (Color online) (a) The effective interaction V(q,oj) 
(double wavy line) between the bosons. The dashed line is 
the bare Bose-Bose interaction 7b, single wavy lines are the 
Bose-Fermi interaction 7bf, and the solid lines are the Fermi 
Green’s function, (b) The Bose self-energies £ n ( g , oj) and 
£ 12 ( 3 , 0 ;). The dotted lines are excitations in and out of the 
BEC. 



Here x(g,a;) is the density-density response function for 
the fermions with momentum q , and frequency oj. The 
corresponding Feynman diagram for V ( 3 , oj) is given in 
Fig. a (a). The momentum dependence reflects the long- 
range of the interaction, as density perturbations prop¬ 
agate in the Fermi gas. Similarly, the frequency depen¬ 
dence of the interaction is due to the fact that it is not in¬ 
stantaneous since these perturbations have a finite speed. 

In the weak coupling BCS limit, k F a F —> 0_ with 
k F the Fermi momentum of the Fermi gas, the density- 
density response function is given by 


X(<bw) 


uf \ 2 Af(eF) 

W 3 [(^) 2 -l] 


(5) 


for frequency/momenta close to the Anderson- 
Bogoliubov sound mode oj = c s q. The velocity is 
c s = v F ^\ + 2 fc F ap/ 7 r/ y/3 [MH23, the density of states 
at the Fermi level e F = k F /2m F is Af(eF) = m F k F / 7 r 2 , 
and z>f = k F /m. In the BEC regime, k F a F —> 0+, 
the Fermi gas becomes a BEC consisting of diatomic 
molecules (dimers) with mass 2m F and density n F /2, 
where n F = k F /3n 2 is the total density of the fermions. 
The density-density response function is then from 
Bogoliubov theory given by [2P1 


x(tw) 


n F g 2 

4 m F (w 2 — OJ q ) 


f *’f \ 2 A/^cf) 


( 6 ) 


Here, oj 2 = q 2 (Am F )~ 1 (q 2 /4m F + 0.67 F 71 F /2) is the Bo¬ 
goliubov spectrum of the dimer BEC, where we have used 
that the scattering length between the dimers is 0 . 6 ap in 
the BEC limit [3Dj. The second equality in © follows 
from the fact that oj q ~ c s q for small momenta, where 
c s = ^/ 0 . 6 aF?TF 7 r/ 2 mp is the Bogoliubov sound speed. 

In general, the density-density correlation function of 
the Fermi gas has a pole at w = c s q in the whole 
BCS-BEC crossover, where c s is the velocity of sound 
for a given scattering length —00 < of < 00 . It fol¬ 
lows from Q that the induced interaction between the 


bosons has the same pole structure: it is attractive for 
oj < c s q , repulsive for oj > c s q , and it diverges when 
oj = c s q. In addition, it has a non-zero imaginary part 
for frequency/momenta inside the quasiparticle contin¬ 
uum of the Fermi gas. It also follows from ©-© that 
the strength k of the induced interaction scales as 

K = Tbf-^ ( £ f)^t (?) 

r^ 

°s 

which should be compared with the strength 7b of the 
direct Bose-Bose interaction. 

We now examine the effects of the induced interaction 
on the excitation spectrum of the Bose gas. To this end, 
we need to calculate the density-density response func¬ 
tion of the Fermi gas in the whole BCS-BEC regime. The 
density response function y(l, 2 ) is defined as a measure 
for how much the density of the Fermi gas changes at 
point (and time) 1 when a potential perturbation S(f> is 
applied at point 2 : 


X(1.2) 
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We apply a QRPA for calculating the Fourier trans¬ 
form of x(l, 2 ) in the superfluid state [28l I'-TITSTTTI . This 
is the simplest microscopic scheme which recovers the 
Anderson-Bogoliubov mode in the BCS regime, and the 
Bogoliubov mode in the BEC regime. It yields a response 
function of the form 


xOtw) = 


Xo(g,^) 

1 — T F L(q, uj) ’ 


(9) 


where xoidi^) is a four-dimensional vector giving re¬ 
sponse due to quasiparticle excitations in the superfluid, 
and A is a 4 x 4 matrix describing the couplings of the 
densities and the order parameter field. The collective 
modes manifest themselves as poles of the density re¬ 
sponse x(< 7 , oj), i.e. as the zeroes of the determinant 


det [1 — T F L(q, w)] = 0. (10) 


The input parameters needed for the QRPA are the 
chemical potential /i and the pairing gap A of the Fermi 
superfluid, which are obtained self-consistently from BCS 
theory. We have for convergence added a small imaginary 
part ir) = ilO -3 e F to the frequencies, and checked that 
the final numerical results do not depend on 77 , as long 
as rj cf- The details of this QRPA calculation can be 
found for example in Refs. and [ SI] , 

Figure [2] shows the speed of the Anderson-Bogoliubov 
mode as a function of 1 /fcpa, determined by finding the 
frequency oj at which the imaginary part of x(?! w ) is 
maximal for a given momentum q. The value of the 
momentum q needs to be chosen small enough so that 
it probes the linear part of the collective mode branch. 
The speed of sound is then the slope c s = oj/q. In the 
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FIG. 2. (Color online) Speed of sound c s in the two- 
component Fermi superfluid as calculated from the pole of the 
density response function x(q, lo). It approaches c s = vy/\/3 
in the BCS limit (the constant green dotted line), and the 
Bogoliubov result c s = VY^/h-ea?/2>n in the BEC limit (blue 
dash-dotted curve). 
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FIG. 3. (Color online) The blue and red regions show 
Rey(g,cj) for /cyo-y = —1. The black solid line is the weak 
coupling Anderson-Bogoliubov mode, and the quasiparticle 
continuum for to > 2A is indicated by a dashed region. The 
green dashed line is the Bogoliubov spectrum e k of the atomic 
BEC, and the red solid line is the Bcliaev spectrum Ek for the 
coupled Bose-Fermi mixture. The damping 5k of the Bcliaev 
excitations is shown as a red dash-dotted line. 


BCS limit, the speed of sound approaches the weakly in¬ 
teracting limit vp /\/3. The numerically calculated speed 
of sound deviates slightly from this in the very weakly 
interacting regime, due to the difficulty of determining 
the slope when the pairing gap is very small. Our nu¬ 
merics reproduce to an excellent accuracy the speed of 
sound results in Ref. [35]. Note that this theory is of 
course not quantitatively correct in the whole BCS-BEC 
crossover. For instance, the speed of sound approaches 
c s = VFy/kpap/3TT in the BEC limit, see Fig. |2l This 
corresponds to a molecular BEC with a scattering length 
2of, instead of the correct value 0.6cif. We emphasize 
however, that the effects discussed below are completely 
general and do not depend on which approximate theory 
we apply to describe the strongly correlated system. 


Figure [3] depicts the real part of the calculated 
density-density response for kpa = —1. At low fre- 
quency/momenta, we clearly see a sharp Anderson- 
Bogoliubov mode where Rex(q, to) changes sign. The dis¬ 
persion of this mode is close to the weak-coupling result 
c s = vp / \/3. For higher momenta, the dispersion curves 
downwards when it approaches the quasiparticle contin¬ 
uum starting at energies above 2A, with A ~ 0.21 ep. 

Once x(<?, w) is calculated, we use Beliaev theory [36] to 
describe the effects of the resulting induced Bose-Bose in¬ 
teraction on the excitation spectrum of the atomic BEC. 
The single particle propagator G(k, w) for the BEC is a 
2x2 matrix, and the Dyson equation reads 

G(k, to) = G 0 (k, to) + G 0 (k, w)£( k, w)G(k, o). (11) 


The bare propagator is 


G 0 (k,w) 


G 0 (k,w) 0 

0 G 0 (k, — w) J 


and the self-energy is 


£(k,w) 


£n(k,w) Si2(k,w) 
S 2 i(k,w) £11 (k,— to) 


( 12 ) 


(13) 


where we have used the inversion symmetry k o —k. 
The effects of interactions are included via the ’’Hartree- 
Fock” self-energies illustrated in Fig. [I] (b), and given 
by £n (k, to) = £n(k, —to)* = n o R(0,0) + n 0 V{k,o) 
and £12 (k, to) = £ 2 i(k,w) = noR(k,u;). Solving these 
equations for G(k, to) yields the Green’s functions for the 
diagonal elements 


G(k, w) 


o; + e fc +n 0 E(k,a;) 
u) 2 -E(k,oj) 2 


(14) 


where E(k,io) = ej" + 2efc?iBE(fc, w). The off- 
diagonal elements are Gi 2 (k, lo) = G 2 i(k, 10 ) = 

—riBE(k,a;)/ [ w 2 — E(k,io) 2 ], where ub is the density of 
the BEC. The theory satisfies the Hugenholtz-Pines re¬ 
lation for the chemical potential /i = £h(0)-£ 12 (0) = 
?4bE(0,0). These interacting Green’s functions describe 
excitations with energy dispersion Ek given by solving 


E k =ReE(k,u = E k ). (15) 


In the absence of the induced interaction, this results in 
the usual Bogoliubov dispersion e k = \/e k + 2n,BTB£k- 
However, due to the momentum and frequency depen¬ 
dence of V(k,oj), (15) is implicit and needs to be solved 
numerically. The equation also yields damping of the 
excitations given by 5 k = lmE(k, E k ). 

Figures i - E show the dispersion E k obtained 
from in the BCS (kpap = — 1), unitarity 1/kpap = 
0, and BEC (kpap = 1) regimes of the Fermi gas re¬ 
spectively. The calculations are performed using param¬ 
eters corresponding to: densities np = Ub = 10 13 cm” 3 , 
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FIG. 4. (Color online) Similar plot as Fig. [3] but for a unitary 
Fermi gas with kpciF = 00. Here A « 0.69 Ef- 



FIG. 5. (Color online) As Figs. [3] and [4] but for the Fermi 
superfluid in the BEC regime with kF<iF = 1- Here A « 
1.35 e F . 


scattering lengths as = obf = 400 ag, and inspired by 
the superfluid Bose-Fermi mixture experiment |24| we use 
the masses of 6 Li and 7 Li atoms. From |7|, this yields 
k = 47 raindmB _ 1 Up/Cg for the strength of the induced in¬ 
teraction with the effective scattering length a; n( j ~ 70 ao- 
Consider first the BCS regime with kpap = — 1 shown 
in Fig. [3j Comparing the Bogoliubov spectrum e*, for 
the atomic BEC decoupled from the Fermi gas with the 
Beliaev spectrum Ek for the coupled Bose-Fermi mixture 
obtained from (15), we see that coupling to the Anderson- 
Bogoliubov mode results in an avoided crossing. Since 
we are neglecting backaction effects on the Fermi gas, 
this avoided crossing becomes a discontinuous jump in 
the bosonic excitation frequency. We expect this predic¬ 
tion to be qualitatively correct, except very close to the 
avoided crossing, since the induced interaction diverges 
when the two excitation frequencies are equal, making 
the corresponding avoided crossing sharp. Figure [3] also 
shows that the excitations of the BEC become damped 
when their energy is inside the quasiparticle continuum 
of the Fermi gas. This reflects that the excitation dissi- 
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FIG. 6. (Color online) Top: correction ( Ek — £k)/t f to the 
BEC dispersion due to the induced interaction as a function of 
Fermi-Fermi scattering length fcFtiF and momentum k. Bot¬ 
tom: decay 8k/t f of the excitations in the BEC. 


pates energy by exciting quasiparticles in the superfluid 
Fermi gas. 

Figure [4] depicts the spectrum Ek when the Fermi gas 
is in the unitarity regime with 1/kpap = 0. We again 
see that there is an avoided crossing, evidenced by a 
jump in the Beliaev dispersion Ek , when the Bogoliubov 
mode approaches the collective mode of the Fermi gas. In 
fact, the resulting energy shift is larger than in the BCS 
case, since the spectral weight of the collective mode is 
larger in the unitarity regime. The bosonic excitations 
are again damped for energies u> > 2A. The small resid¬ 
ual damping near the avoided crossing reflects however 
the small imaginary part irj that we have built into the 
Fermi theory to obtain convergence. In the limit rj —> 0, 
the bosonic excitations are undamped outside the quasi¬ 
particle continuum, even at the avoided crossing since it 
corresponds to the coupling of two undamped excitations. 

Finally, Fig. [5] shows the dispersion Ek in the BEC 
regime of the Fermi gas with kpap = 1. The avoided 
crossing feature and the energy shift in Ek is now even 
more pronounced due to a smaller sound velocity of the 
Fermi gas, which approaches the Bogoliubov sound speed 
of a dimer BEC, thereby making n larger as can be seen 
from 0. The quasiparticle continuum of the Fermi gas 
is outside the range of the plot due to the large pairing 
energy in the BEC regime. There is therefore no damping 
of the bosonic modes shown. 

The above results show how the coupling between the 
superfluid bosons and fermions leads to significant ef¬ 
fects on the spectrum of the atomic BEC, which depend 
on the properties of the Fermi gas. In the recent exper¬ 
iment on the superfluid 6 Li and ' Li mixture, the Fermi- 
Fermi scattering length a F could be tuned using a Fesh- 
bacli resonance. We therefore plot in Fig.[ 6 ]the difference 
Ek — Sk between the Beliaev and Bogoliubov excitation 
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spectra as a function of czf, keeping all other parame¬ 
ters as in Figs. 3][5 We also plot the damping of the 
mode. Two effects are apparent. First, since the sound 
velocity in the Fermi gas depends on ap, the momentum 
where the bosonic mode exhibits the avoided crossing 
depends on of- Also, the induced interaction in gen¬ 
eral decreases/increases Ek for energies below/above the 
avoided crossing as expected. The magnitude of the en¬ 
ergy shift increases towards the BEC regime since spec¬ 
tral weight of the collective mode in the Fermi super¬ 
fluid increases. Second, Fig.[ 6 ]clearly shows the damping 
caused by the coupling to the quasiparticle excitations 
of the superfluid Fermi gas. This quasiparticle contin¬ 
uum moves to higher momenta as the system approaches 
the BEC limit and the pairing gap increases. The resid¬ 
ual damping below the quasiparticle continuum shown in 
Fig. [ 6 ] is, as explained above, a result of using a non-zero 
77 in the numerics, and it vanished for 77 —>• 0. This il¬ 
lustrates how the collective and single particle spectrum 
of the strongly correlated Fermi gas can be mapped out 
by measuring its effects on the excitations in the BEC. 
We note that the effects can be increased significantly 
by increasing a bf, since k oc a| F . In addition to vary¬ 
ing of and obf, one can also vary as which will increase 
even further the ways one can probe the excitations in 
this Bose-Fermi mixture. The excitations of a BEC have 
already been measured using Bragg spectroscopy EMU- 

In conclusion, we examined a mixture of a BEC and a 
superfluid Fermi gas using Beliaev theory for the bosons 
combined with quasiparticle random phase approxima¬ 
tion for the fermions. The fermions were shown to me¬ 
diate a frequency/momentum dependent interaction be¬ 
tween the bosons, which leads to two qualitatively new ef¬ 
fects. First, the induced interaction diverges at the sound 
mode of the Fermi gas which results in a sharp avoided 
crossing feature in the excitation spectrum of the BEC. 
Second, the excitation of quasiparticles in the Fermi gas 
leads to a damping of the excitations of the BEC. By 
varying the densities and scattering lengths of the sys¬ 
tem, these effects can be used to systematically probe 
the properties of the Fermi gas in the strongly correlated 
BCS-BEC cross-over. Our work may be extended in a 
number of directions: It would be interesting to include 
the backaction of the bosons on the superfluid Fermi gas 
to obtain a detailed description of the avoided crossing 
of the sound modes. Trapping effects can be included 
using a local density approximation, which has proven 
to work well when considering short wavelength Bragg 
scattering [35]. Finally, the theory can be extended to 
finite temperatures, which would result in a damping of 
BEC excitations for all momenta due to the presence of 
thermally excited quasiparticles in the Fermi gas. 
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